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ABSTRACT 



Aims. The detection of close and faint companions is an essential step in many astrophysical fields, including the search for planetary 

companions. A new method called "phase closure nulling" has been proposed for the detection of such faint and close companions 

based on interferometric observations when the system visibility amplitude is close to zero due to the large diameter of the primary 

star We aim at demonstrating this method by analyzing observations obtained on the spectroscopic binary HD 59717. 

Methods. Using the AMBER/VLTI instrument in the K-band with ~ 1500 spectral resolution, we record the spectrally dispersed 

closures phases of the SB 1 binary HD 597 17 with a three-baseline combination adequate for applying phase closure methods. After a 

careful data reduction, we fit the primary diameter, the binary flux ratio, and the separation using the phase closure data. 

Results. We detect the 5-mag fainter companion of HD 59717 at a distance of 4 stellar radii from the primary. We determine the 

diameter of the primary, infer the secondary's spectral type and determine the masses and sizes of the stars in the binary system. This 

is one of the highest contrasts detected by interferometry between a companion and its parent star 
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1. Introduction 



stars: individual: HD 5917 - binaries: close - binaries: spectroscopic - techniques: 



Since the first firm discoveries of planetary-mass companions 
to stars, there has been a renewed interest in all the observa- 
tional techniques that can provide reliable estimates of binary 
. masses and distances. Beside the mass, the most desirable in- 

■ formation and the most difficult to measure is the spectrum 
of the companion. Imagi ng the close environment of stars is a 

. very activ e reseai'ch field (iMarois et al . 2008: Kalas et al.ll2008t 
' iLagrange et al.. ,20091) . and is the main m otivation for several 
'. ambitious instrume nts proposed recently (ICockell et al.l 20091 
iLawson et al.'"2008'), based either on fringe nulling ( Bracewelll 
1978; Woolf& Angel 1998^ or "extreme" adaptive optics and 
, coronography (see review bv Beuzit et al.ll2007l) . These tech- 

■ niques are complementary to interferometric imaging since they 
I are bUnd to regions within a few Airy disks from the central star, 

whereas the interferometric observables acquired by spatially- 
filtered interferometers largely come from within this Airy disk. 

A large fraction of the observational effort of current op- 
tical long-baseline interferometers is devoted to the study of 
the close environment of stars. By using spatial filtering tech- 
niques, the interferometers have gained in precision and stabil- 
ity. Using more than two apertures simultaneously, together with 
a large spectral coverage, near infrared beam combiner like the 
AMBER/VLTI instrument d Petrov et aT]|2007l) . can provide pre- 
cise measurements of these close environments, of less than a 
few hundred stellar radii. We have recently proposed a new inter- 



ferometric technique called "phase closure nulling' (hereafter 
PCN), to detect and to perform the spectroscopy of faint com- 
panions of stars (Chelli etal. 2009a). This technique is based 
on the modeling of spectrally dispersed phase closure measure- 
ments of the multiple system around visibility zeros of the pri- 
mary. In these regions, there is always a spatial frequency inter- 
val within which the phase closure signature of the companion 
is larger than any systematic error and is thus measurable. 

We illustrate the proof of concept of PCN with a simple 
observational case, the bright single-lined spectroscopic binary 
HD 597 17. Although the observations were not initially intended 
for this purpose, and are thus very incomplete, we can neverthe- 
less derive the characteristics of the system in terms of stellar di- 
ameter, separation and flux ratio. The observations and the data 
reduction are presented in Sect.|2] Section[3]describes the deriva- 
tion of the binary parameters, from visibility and phase closure 
measurements. The results are discussed in Sect.|4l 

2. Observations and data reduction 

2.1. Observations 

The observations are part of a series of on-sky tests performed 
to assess the stability of the AMBER instrument and the corre- 
sponding accuracy of absolute calibration. They were performed 
after a change of the infrared detector and the temporary re- 
moval, by our team, of a set of polarisers that was inducing 



* Based on observations collected at the European Southern 
Observatory, Paranal, Chile, within the commissioning programme 
60.A-9054(A). 



' Not to be confused with the "Three-Tele scope Closur e -Phase 
Nulling Interferometer Concept" proposed by iDanchi et al.l l l2006h 
which deals with the phase closure properties of nulling interferome- 
ters. 
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Table 1. Summary of observation log. Uniform disk diameters for calibrators were obtained with the use of SearchCal 
(iBonneau et al.ll2006l) . 



Name 


Date 


Time of observation 


Sp. Type 


Mag. K 


^UD 


Fringe 


Spectral coverage 


DIT 












(mas) 


Li aLKCi 


(nm) 


(ms) 


HD 59717 


2008-02-13 


03:03:07, 05:40:59 


K5III 


-0.41 




yes 


[1925-2275] 


200 


HD 56478 


2008-02-13 


02:26:23, 04:26:38 


KOIII 


4.48 


0.618 ±0.043 


yes 


[1925-2275] 


200 


HD 37350 


2008-02-13 


05:04:59 


F6Ia 


2.09 


0.643 ± 0.026 


yes 


[1925-2275] 


200 


HD 59717 


2008-02-14 


03:02:50,04:51:56, 06:11:50 


K5III 


-0.41 




no 


[2110-2190] 


50 


HD 70555 


2008-02-14 


03:27:16, 05:18:42, 06:40:03 


KlIII 


1.53 


2.54 ± 0.037 


no 


[2110-2190] 


50 


HD 45669 


2008-02-14 


03:49:44, 04:26:21, 05:48:39 


K5III 


1.88 


2.19 ±0.053 


no 
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Fig. 1. Transfer function computed from the 6 calibrator obser- 
vations carried out on 14 February 2008 as a function of time 
for the 3 baselines (shortest baseline at the top). The numbers in 
parenthesis are the relative visibility dispersion. 



Fabry-Perot fring e beating, leading t o instabilities in the instru- 
mental visibility (IMalbet et al.ll2008h . As a consequence of this 
removal, the fringe contrast decreased, but the stability improved 
drastically. To characterize the stability, we observed a series of 
calibrator stars of various magnitudes, at different locations in 
the sky, taken from the ESO list of VLTI calibrators. One of these 
stars, HD 59717, is a triple system located at 56.3 pc, consisting 
of a K5III single-lined spectroscopic binary and a G5V compan- 
ion at 22'.'3 distance (MSG catalog, Tokovinin 1997). The com- 
panion was not in the field of view of the AMBER observations, 
and we will refer hereafter to the SBl close binary as HD 59717. 

We used three 1 .8m auxiliary telescopes (ATs) on VLTI sta- 
tions HO, DO and AO. In this configuration, the 3 telescopes are 
aligned in a very rough East-West (71°) configuration and the 
ground baseline lengths are AO-DO: 32 m, DO-HO: 64 m, AO-HO: 
96 m. Due to its large angular size, HD 59717 is well resolved 
in the HO-DO-AO configuration, the longest baseline crossing 
by supersynthesis effect the first zero of the visibility curve and 
providing information in the first and the second lobes. 

Our set of data consists of three observations of HD 59717 
performed on 14 February 2008, using the spectral window 
2110-2190nm and a detector integration time (DIT) of 50 ms. 
Each observation was bracketed with those of the calibrators 
HD 70555 and HD 45669. Additionally, we make use of a sin- 
gle visibility measurement taken from an observation performed 
the pr evious night with FINITO, the VLTI fringe tracking fa- 
cility (iGai et a l. 2004). This observation covers the wavelength 
range 1925-2275 nm, used a DIT of 200 ms, and was bracketed 
with observations of the calibrators HD 56478 and HD 37350. 



2.2. Data reduction 

Prior to the data reduction, we performed an accurate wavelength 
calibration thanks to the identification of atmospheric features in 
the spectrum obtained from the single observation with FINITO. 
We fit the position of the atmospheric features with a second 
degree polynomial, providing a precision of a few nanometers. 
Then, the data were redu ced with the new A MBER data reduc- 
tion method presented in lChelli et"an (l2009bl) . 

Figure [1] shows the instrumental transfer function averaged 
in wavelength, obtained for the 3 baselines as a function of time. 
This transfer function results from the ratio of the measured visi- 
bility and the intrinsic visibility of the calibrator. The number be- 
tween parenthesis is the relative visibility dispersion, computed 
for each baseline, as the ratio between the dispersion of the vis- 
ibilities along the night and their mean. This is the dominant 
source error due to a residual atmospheric jitter during the frame 
integration time. Other sources of error are negligible in this re- 
spect. The intrinsic visibilities have be en derived from diam eters 
estimated with the SearchGal software dBonneau et al .1200 6') and 
their associated relative error is much smaller, less than 0.8% for 
the largest calibrator on 100m baseline. The high spectral res- 
olution (1500) prevents any spectral mismatch between source 
and calibrators. Except for the single visibility data in the sec- 
ond lobe obtained with FINITO (see below), the two magnitude 
difference between the source and the calibrators does not matter 
as, in the absence of a fringe-tracking system, the jitter is flux- 
independent. In the spectral range 2110-2190nm, the transfer 
function is a linear function of the wavelength. We use this a 
priori information to reduce the high frequency noise on the cal- 
ibrator visibilities by Fourier filtering. 

The visibiUty data of HD59717 obtained with the assistance 
of the fringe tracker (see Table [TJ are affected by a residual, 
flux-dependent, jitter that will be different on our source and its 
calibrators due to their magnitude difference. In principle, we 
can apply a jitter correction to the data. Unfortunately, the jit- 
ter residuals from the fringe tracking system were not included 
at the time in the instrumental data. Only one visibility of high 
quality, measured in the second lobe, could be calibrated with the 
aid of partially redundant visibilities obtained during the second 
night. 

The peculiarity of AMBER is to use an internal calibration 
called P2VM to estimate the coherent fluxes. Hence, the phase 
closure on a source contains the imprints of the P2VM. Since, 
during the building process, the three beams are measured in- 
dependently, the internal phase closure of the P2VM is not by 
construction constrained to be zero. This is why a phase closure 
measurement with AMBER must be calibrated by substracting 
that of a calibrator The three (spectrally dispersed) phase clo- 
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Fig. 2. Visibility of HD 59717 as a function of the spatial fre- 
quency. Errorbars illustrate the calibration errors resulting from 
the visibility dispersion shown in Fig.[T] The measurements ob- 
tained with the fringe tracker are plotted in grey in the second 
lobe. The full curve corresponds to a single uniform disk model 
of diameter 6.436 mas. The insert shows the second lobe mea- 
surements and three binary models with separation 30 mas and 
flux ratio: 5 x 10"^ (continuous line), 1 x 10"^ (dashed line) , 
2 X 10^2 (dotted line). 



sures measured just before, during and after the minimum of vis- 
ibility have been corrected by the phase closure of a calibrator 
star close in time. The main error, due to the calibration process, 
is about 0.05 radian from the internal dispersion of the calibrator 
phase closure measurements along the night. To reduce the phase 
closure dispersion, especially at minimum visibility, we average 
the bispectra (from which the phase closure is computed), ini- 
tially sampled with a frequency interval of 0.06 arcsec" ' , whithin 
a frequency interval of 1 arcsec 



3. Results 

The calibrated visibilities of HD 59717 are shown in Fig.|2]as a 
function of the spatial frequency, defined as the ratio between the 
baseline and the wavelength. The measurements are reasonably 
uniformly spread over the working frequency range, with a set 
of measurements before, partially during, and after the minimum 
visibility. The full curve represents the best fit (excluding the 
points around minimum of visibility, see Sect. 13. lb with a single 
uniform disk model of diameter 6.436 mas. This is in excellent 
agreement with the diameter of 6.44 ± 0.06 mas arising from the 
location of the minimum of visibility around 189 ±0.02 arcsec"^. 

The phase closures are shown in Fig. [3] The transition be- 
tween the and n values at the zero visibility crossing is smooth, 
markedly different from the expected abrupt step arising from 
a cent rally symmetric flux distribution only (see IChelli et"al] 
l2009ah . 

This departure from a centrally-symmetric object is expected 
in a spectroscopic binary, and we can use these observations to 
retrieve the geometric characteristics of the system. To do this, 
we fit independently the visibility and the closure phase data 
with a simple model formed by an extended uniform disk and 
a point source. We did not try to combine visibility and closure 
phase data because the weight to give to each set of data is quite 
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Fig. 3. Phase closure of HD 59717 as a function of the largest 
spatial frequency of the closure triangle. The thin curve cor- 
responds to the best fit for a single uniform disk model with 
diameter 6.451 mas. The full curve corresponds to the best fit 
with a double system and parameters: primary stellar diameter 
6.55 mas, secondary projected distance - 11 .2 mas and flux ratio 
0.017. 



uncertain. Our observations being basically east-west, we can- 
not derive the position angle of the system, but only the pro- 
jected distance along the direction of observation. As we have 
used the spectral dimension to increase the frequency coverage, 
we assume that the object is achromatic. Hence our model is de- 
scribed by 3 parameters: the stellar diameter, the distance of the 
companion and the flux ratio. 



3.1. Visibility data 

To estimate the system parameters, we minimize a standard 
defined as the distance between the visibilities of the model and 
those of the data, weighted by the errors shown in Fig. [T] We 
exclude from the fit the set of points around minimum visibility 
as they may easily be biased. Indeed, at this location, the bias on 
the squared amplitude of the coherent flux is twenty times larger 
than the useful (debiased) signal. Given the relative error of 0.07 
on our longest baseline, an imprecision of 0.07/ V20 =1.5 10"^ 
would produce a bias "error" equivalent to the statistical error. 
Since we are not guaranteed such a precision in the bias removal, 
it is better to exclude this set of points, which otherwise would 
constrain nearly by itself the output of the fit. This restriction 
does not apply if a single stellar diameter is deduced from the 
position of the minima of visibility, since the location of the min- 
ima are much less biased. 

A rapid study of the problem shows that the has min- 
ima as a function of the distance with a peudo-period p x /^^.' , 
where /~/^ is the frequency of minimum visibility. Hence to pro- 
duce comprehensible output, we vary the separation from to 
0.1 arcsec, and for each separation, we perform a fit with two 
parameters, the stellar diameter and the flux ratio. The minimum 
reduced;^^^ (hereafter ;i'^), the best stellar diameters and flux ratio 
are displayed in Fig.|4]as a function of the separation. 

The;\f^ presents shallow minima between 1.3 (single source) 
and 0.9 (double source) at a set of roughly regularly separations. 
However, given that below 1 % precision, the calibration biases 
on the visibility become dominant, the presence of these minima 
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Fig. 4. Best model parameters from visibility fit with an extended 
uniform disk and a point source (see Sect. l3.1l for details). From 
top to bottom as a function of the separation: minimum best 
flux ratio and best stellar diameter. 



alone does not prove the binary nature of HD 59717. In addi- 
tion, their location must be taken with extreme caution. If we 
except the singular point around 3 mas inside the stellar disc, 
the range of possible stellar diameters for all the separations is 
6.45 + 0.03 mas. This diameter range agrees with the output of 
the fit with a single uniform disk model of 6.436 mas (see Fig.|2]). 
The possible flux ratio ranges from 10"^ to 6 x lO""' for all sep- 
arations and, the smaller the separation the higher the flux ratio. 



3.2. Phase closure data 

The best fit of the phase closure systems with a single uniform 
disk model provides a diameter of 6.451 mas in excellent agree- 
ment with the value of 6.436 mas derived from visibility data. 
The minimum of a two-component model, the best stellar di- 
ameters and flux ratio are displayed in Fig.|5]as a function of the 
separation. This minimum varies from 120 (single source) to 
2.3 (double source). 

The stellar diameter ranges between 6.2 and 6.7 mas. It is 
not well constrained because we do not have the full transition, 
but only pieces of it. The flux ratio varies between 6 x 10"^ and 
2.7 X 10"^, which corresponds to a 5^y^g mag difference between 
the primary and the secondary. The slope of the phase closure is 
extremely sensitive on the flux ratio. To illustrate this effect, we 
plotted in Fig. |6]the phase closure together with the best fit for 
various flux ratios. One sees that all flux ratios < 5 x 10"^ cannot 
reproduce the slope of the transition. 

The;^^^ (Fig.|5]l exhibits regularly spaced and very deep min- 
ima as a function of the separation. The minimum flux ratio 
of 5 X 10"^ obtained in the fits renders separations larger than 
30 mas irrelevant, as they would produce oscillations that are not 
seen in the visibility, especially in the second lobe (see insert in 
Fig. |2]i. Hence, the possible projected separations left are -5.5, 
-1-8.5, -11, -(-14 and -17 mas, with an error of ~ 1 mas, from the 
width at mid-height of the;^'^ dips. 

In summary, visibility data put strong constraints on the stel- 
lar diameter, 6.45 ±0.03 mas. Phase closure data do not constrain 



Fig. 5. Best model parameters from the phase closure fit with an 
extended uniform disk and a point source (see Sect. l3.2l for de- 
tails). From top to bottom as a function of the separation: mini- 
mum jt'^, best flux ratio and best stellar diameter The curves rep- 
resent the best fit with a restricted stellar diameter in the range 
6.42-6.48 mas imposed by visibility modelling. 




Fig. 6. Best fit of the phase closure for different flux ratios show- 
ing that the slope of the transition cannot be reproduced with a 
flux ratio smaller than 5 x lO^^*. 



the stellar diameter but support the binary nature of HD 59717 
(given the impossibility of reproducing with a single uniform 
disk model the slope of the phase closure around the minimum 
of visibility), providing a range of acceptable flux ratios and five 
possible separations. 

4. Discussion 

The phase closure transition of HD 59717 is markedly different 
from that of a single stellar disk. Our modeling of this phase clo- 
sure transition is naturally based on the a priori information that 
the system is a binary. We do not consider here the influence of 
large stellar sp ots because K giants are devoid of such features 
(ICohen et all 1999) . and because the larger of solutions inside 
the stellar radius in Fig. |5] make it unlikely. Although we detect 
the faint companion by its effect on the phase closure, the very 
incomplete information we have on its shape around the first null 
(compare our Fig.|3]with Fig. 3 of IChem etal.ll2009al) prevents 
a precise determination of the secondary's position using the in- 
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Table 2. Literature data for HD 59717 (from lJancart et al.ll2005HWilsodll918l:IPerrvman & ESAlll997h . 
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Fig. 7. Loci of the solutions for the equation of mass (black), as a 
function of the mass mi of the primary and m2 of the secondary. 
The regions compatible with the set of 3 separations given by 
our fit are overlayed in gray. The region compatible with our 
independent measurement of the secondary's flux is overlayed 
in light gray. 



terferometric measurements alone. It is thus important to check 
whether our results are in agreement with the already known pa- 
rameters (star diameter, spectroscopic orbit) of HD 597 17. 

Table |2] summarizes the literature data for HD 59717 rel- 
evant to this discussion. HD59717, at a distance of 56.36pc, 
has a spec troscopic peri od of 257'?8 and a reduced mass yU - 
O.164M0 (IWilsonlll918h . The parameters of the photometric or- 
bit of HD 59717, in particular the inclination / = 65% + 3°.3 
and semimajor axis ap - 8 .32 ± 0.32 mas, have been ob- 



tained by Jancart et al.| (|2005 ) from Hipparcos IAD measure- 
ments (Ivan Leeuwen & Evans I998I) . HD 59717 exhibits a small 



amplitude (~ 30millimag in V) intrin sic variability character- 
istic of rota ting ellipsoidal variables (iPerrvman & ESAI Il997t 
IOteroll2008l) . 



4.1. Stellar diameter of the primary 

The stellar diameter of HD 59717 has been estimated by 
ICohen et al.l (Il999l) who derive a limb-darkened (LD) angular 
disk diameter of 6.86 ± 0.1 mas. Using the relation be tween the 
LD and UD diameters of lHanburv Brown et alj (|I974|) and a lin- 
ear li mb-darkening coefficient at 2.2//m of 0.379 (ICIaret et all 
119951) . this translates into an equivalent UD diameter in the K 
band (UDK) of 6.66 + O.I mas. The SearchCal web service 
dBonneau et al ] l2006h gives UDK=6.16 ± 0.4 mas. Our observa- 



tions show that the diameter of the primary is well constrained 
by the visibility data alone to UDK=6.45 + 0.03 mas. 

4.2. Nature of the companion 

In stellar systems for which all the orbital and projection pa- 
rameters are measured, the knowledge of the mass mi of the 
main component is sufficient to deduce the mass, hence the 
spectral type, magnitude, and true orbit of the unseen compan- 
ion star Given mi, the reduced mass and the inclination /, 

one solves the equation of mass: ^"^"^^^'2 = fi for the mass m2 
of the secondary. When the primary is a main-sequence star, 
mi is reli ably estimated using a stand ard mass-luminosity func- 
tion (e.g. jMcCIuskev & Kondoll972h . Unfortunately, there is no 
such thing as a unique mass-luminosity function for giant stars 
(as is the case for HD 59717), so mi is unknown and individ- 
ual masses in a binary system with a giant component cannot be 
derived reliably by this method. 

Conversely, if the companion of the giant star is a main- 
sequence star, and if its luminosity can be measured, then m2 
is known and the mass mi of the giant can be measured. In the 
present case, our observations provide two new pieces of infor- 
mation: the luminosity of the secondary via our measurement of 
the flux ratio, and a choice of projections of the separation on 
the E-W direction. Each of these pieces of information can be 
used independently to derive the mass of the giant primary, and 
should provide compatible results. 

4.2.1. Position constraints 

We first combine the constraints due to the measured projected 
separation and the equation of mass. In Fig. [T] we show the loci 
of the solutions for the equation of mass with / = 65°6 ± 3°3, fi = 
0.164Mq in black. At the time of our observations, the position 
angle of the binary was 6 - -76°0 + 1°5, and the projection of 
the separation sp - a(l + '■^), where a - -3.5 ± 0.2 is fixed by 
the orbit geometry. The positive projected separations (+8.5 and 
+ 14 mas) are forbidden by the value of a. The three possible 
(ml, ml) solutions left are for the projected separations -5.5, 
-11 and -17 mas, and are plotted in Fig. |7]in gray. Their width 
corresponds to the uncertainty (+1 mas) on the separation. 

The solutions for the binary masses lie at the intersection be- 
tween these areas and the curve given by the equation of mass. 
The -5.5 mas and -17 mas regions give unrealistic values for 
the mass distribution in the binary: the former since the primary 
would be less massive than the secondary, the latter because 
it gives a too high mass for the primary. This leaves only the 
-11 mas solution, and our mass estimate by this method is then 
a 5M0 primary and a 2.2M0 secondary. 

4.2.2. Pliotometry constraints 

If we use now the magnitude of the secondary derived from our 
observed flux ratio, its spectral type is A0 V-A2V (see Table[3]l. 
Then mi lies between 2.0 Mq and 2.5 Mq (IMcCIuskev & Kondol 
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Table 3. Known (bold text), measured (italic) and tabulated parameters for the A and B components of the HD 59717 binary. 



Star 


mK 


MK 


AMk 


MV 


AMv 


Sp.Type 


Mass 


Radius 




log g/g0 


logp/po 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




(Mo) 


(mas) 


(Ro) 






A 


-0.41 


-4.16 




-0.55 




K5III 


5 ±0.5 


3.23±0.015 


39 


-2.5 


-4.7 


B 


5.14 


1.39 


5.5 


1.8 


2.4 


A4V 


2.0 


0.15 


1.8 






B 


4.89 


1.14 


5.3 


1.3 


1.9 


A2V 


2.2 


0.16 


2.0 






B 


4.59 


0.84 


5.0 


0.65 


1.2 


AOV 


2.5 


0.2 


2.4 







[T971 . In Figure |7] we have also reproduced this secondary's 
mass range as a horizontal area shaded in light gray. It inter- 
sects the curve of the equation of mass at 4.5 < m\ < 5.5 M0, 
and thus independently gives the same result as our positional 
measurement. 

Table[3]collects all the parameters on the primary that can be 
deduced from our size and mass measurements. Three possible 
identifications are proposed for the B companion, allowing for 
the uncertainty on the flux ratio. The sectral types and radii of 
the companion are fromiSchmidt-Kaler (1982), mas s es fro m the 
mass-luminosity relation of iMcCluskev & Kondol (|I972|) . We 
note that the radius and the mass of HD 59717 are larger than 
th e values {Rksiii - 25 R© and Mkhii = 1-2 M©) published 
bv lSchmidt-Kalerl ([T982h . However, the masses we obtain are in 
agreement with the v alues (Maa - 5 M ©, M^z, > 1 .9 M©) quoted 
in the MSC catalog (iTokovininI 1 997h for the Aab component of 
HD 59717. 

4.3. The photometric variability of HD 5971 7 

The separation of the components, their size, and the inclination 
of the orbit make it impossible for the secondary to be eclipsed 
by the primary. With a 5 Mq primary and a 2.2 M© companion, 
the Roche lobe Rl radius of the primary is Ri ki 0.45 x d, where 
d is the true distance of the two stars. Due to the excentricity 
e - 0.17 of the companion, this distance varies from 23 mas 
(periastron) to 32 mas (apoastron). With its 3.23 mas radius, the 
primary thus occupies between 22% and 32% of its Roche lobe, 
and is slightly ellipsoidal. This is accordance with its classifica- 
tion as a rotating ellipsoidal variable. 

If the photosphere of the star was of uniform brightness and 
shaped by the Roche potential, the maximum prolateness of the 
star would occur at periastron and he a/b - 1 .0022. At the time 
of our observations, the deviation of the projected shape of the 
star from a circle would be ~ 6 x lO"**, in accordance with our 
hypothesis that the primary can be modeled as a uniform disk. 
However, with such a small prolateness, the change in the pro- 
jected surface of the star with time would be < 10"^, which can- 
not account for the ~ 30 millimag photometric variability, which 
requires a/b ~ 1.02. The variability of HD 59717 could be due to 
a difference of surface brightness between the hemisphere of the 
star facing the companion and the other hemisphere. Although 
this effect should, to a first approximation, translate as a dis- 
placement of the photocenter of the star and would not change 
the phase closure values, it will be investigated more closely in 
a future study aimed at confirming at another epoch, with ded- 
icated observations, the detection of HD 59717 companion by 
phase closure nulling. 

5. Conclusion 

We have shown that, by fitting a simple model on the spectrally- 
dispersed phase closure measurements made by AMBER on the 



bright SBl binary HD 59717, observed near the first zero of 
the visibility distribution, we can resolve the close (~ 4 stel- 
lar radii at the time of observation) pair with a five magnitude 
difference in brightness. This result, obtained with only 15 min- 
utes of on-sky integration and 1.8 m apertures, is two magni- 
tudes fainter than reported detection limits (applying Eq. A. 2 
of Tokovinin et al. 2006, for p = 1 1 mas) and is an illustra- 
tion of t he potential of the p hase closure nulling method pro- 
posed bv Chelli et alJ(l2009ah . We have detailed the specific data 
processing involved, derived the spectral type of the compan- 
ion and measured the individual masses and sizes of the system 
(TableO. 

The PCN method employed here is of interest for the charac- 
terisation of all faint point-like sources in the immediate vicin- 
ity of stars, if the flux from the primary is sufficient and with 
the quite restrictive condition that the primary can be fully re- 
solved by the interferometer. This technique potentially gives ac- 
cess to masses, eve n spectra, of brown dwarfs and exoplanets 
(ICheUi et alJl2009ah . In this paper we have used published astro- 
metric data, since our observations were of incomplete coverage 
both in spatial frequency and position angle, but dedicated PCN 
observations taken at different epochs should suffice to charac- 
terize companions and fit their orbit. For exoplanets detected by 
radial velocity techniques, a single PCN measurement of the sec- 
ondary's position will solve for inclination of the orbit of the 
planet. Moreover, since radial velocity techniques can be con - 
fused by the presence of stellar spots (see .Huelamo et al.ll2008l) . 
PCN can provide the independent detection needed for confir- 
mation of the presence of an exoplanet. 
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